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The device for the Josephson flux qubit (DJFQ) can be considered as a solid state artificial atom 
with multiple energy levels. When a large amplitude harmonic excitation is applied to the system, 
transitions at the energy levels avoided crossings produce visible changes in the qubit population 
over many driven periods that are accompanied by a rich pattern of interference phenomena. We 
present a Floquet treatment of the periodically time-dependent Schrodinger equation of the strongly 
driven qubit beyond the standard two levels approach. For low amplitudes, the average probability 
of a given sign of the persistent current qubit exhibits, as a function of the static flux detuning 
and the driving amplitude, Landau-Zener-Stiickelberg (LZS) interference patterns that evolve into 
complex diamond-like patterns for large amplitudes. In the CEise of highly coherent flux qubits we 
show that the higher order diamonds can not be simply described relying on LZS transitions in each 
avoided crossing considered separately. In addition we propose a new spectroscopic method based 
on starting the system in the flrst excited state instead of in the ground state, which can give further 
information on the energy level spectrum and dynamics in the case of highly coherent flux qubits. 
We compare our numerical results with recent experiments that perform amplitude spectroscopy to 
probe the energy spectrum of the artiflcial atom. 
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I. INTRODUCTION 

In recent years, several types of superconducting qubits 
have been experimentally studiedji"— These systems con- 
sist on mesoscopic Josephson devices and constitute 
promising candidates to be used for the design of qubits 
for quantum computationjir— Indeed, a large effort is de- 
voted to succeed in the coherent manipulation of their 
quantum states in a controlable way. The progress made 
along this line allows to have nowadays Josephson circuits 
with small dissipation and large decoherence times.— "^i^"— 

In this work we will focus on the device for the Joseph- 
son flux qubit (DJFQ), which consists of a SQUID loop 
with three Josephson junctions operated at or near a 
magnetic flux of half quantum. -i2r^. When cooling down 
to millikelvin temperatures this device exhibits quantized 
levels whose energies can be tuned by a control param- 
eter such as an exernal magnetic field. This artificial 
atom-like behavior has motivated several studies based 
on the analysis of the level spectrum and its dynamics 
beyond the simplified two-level approach. As an exam- 
ple it has been shown that, after the inclusion of higher 
energy levels the DJFQ exhibits quantum signatures of 
classical chaoSf2i^ and a recent studyiii focused on the 
calculation of the intrinsic leakage ( i.e. transitions from 
the allowed qubit states to higher excited levels of the sys- 
tem) has shown that for very strong resonant harmonic 
pulses the two-level approximation breaks down. 

What is more important, several recent experiments 
driving the fiux qubit with a combination of a dc and 
large amplitude harmonic excitations in the magnetic 
fiux have studied the energy level structure through 
Landau-Zener-Stiickelberg transitions.—"— Mach-Zender 
interferometrjii^ii^ and amplitude spectroscopj»i^ have 
been the subject of these recent experimental studies 



of the flux qubit as an artiflcial atom. In par ticular, 
the amplitude spectroscopy experiment of Reflla has re- 
vealed the higher energy level spectrum when increasing 
the microwave amplitude. In this case, the average pop- 
ulation of one state of the DJFQ as a function of the dc 
fiux (fiux detuning) and microwave amplitude exhibits 
diamond-like interference patterns, which display a rich 
structure of multi-photon resonancesJ^ From these inter- 
ference patterns it is in principle possible to reconstruct a 
large fraction of the energy spectrum and methods based 
on two dimensional Fourier transform have been recently 
proposed to this end.^^ In the experiment, the observed 
spectroscopic "diamonds" arise due to combined contri- 
butions of Landau-Zener-Stiickelberg transitions, which 
provide the interference fingerprint of different energy 
level avoided crossings, together with intra-well fast re- 
laxation and short coherence times, which provide con- 
trast in the observed pattern. 

Recent theoretical efforts have been put forward to re- 
construct the experimentally observed interference pat- 
terns by solving the dynamics of the model under strong 
driving. Most of the reported approaches reduce the 
model for the DJFQ to a simplified version which only 
considers the dynamics of the two levels involved in each 
avoided crossing. In this case the well known Landau 
Zener- Stiickelberg theory has been applied, consider- 
ing only the accumulated phase of the two levels dur- 
ing a period of the driving] ^^d^ The beginning of the 
first spectroscopic diamond, that corresponds to the first 
avoided crossing, is accurately reproduced within this ba- 
sic modeli ^^d"^ Additionally, extensions that incorporate 
several levels, based on rate equations along all coupled 
levels have been recently proposed. In this later case 
more than one diamond can be obtained, but the ap- 
proach neglects the phase accumulated in the evolution 
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of several levels and can only be applied when decoher- 
ence and relaxation effects are important. Since this case 
is somewhat near the experiment of ReflTsl. a qualitative 
description of the observed diamond patterns can be ob- 
tained. However, the devices for the flux qubit can have 
larger decoherence times^ii^ than in the case ofRefEl. In 
this case, the effect of fast intra-well relaxation that pro- 
vided contrast in the diamond patterns, and the effect 
of short decoherence times that made possible to con- 
sider only the accumulated phases of two-levels at the 
avoided crossings, will be much weaker. Then the ques- 
tion arises on how the interference patterns of strongly 
driven D JFQ can be analyzed in the highly coherent case, 
and how much information on the energy spectrum can 
be extracted in this situation. 

The purpose of this work is to solve the dynamics of the 
DJFQ under strong driving rf pulses considering the full 
hamiltonian of the system. We perform a first-principles 
calculation taking as input only two parameters from 
the experimental device: the ratio of the Josephson and 
charging energies, Ej jEc, and the asymmetry factor a of 
the Josephson energy of one of the junctions with respect 
to the others. Furthermore, our approach focus on the 
behavior of highly coherent DJFQ when driven within 
time scales smaller than the dephasing time. Therefore 
the interaction with the environment is neglected and we 
solve the time dependent Schrodinger equation consider- 
ing the DJFQ as a closed system. As we will show, even 
when relaxation and dephasing are neglected, our results 
reproduce several of the qualitative features of the exper- 
iment of Reffisl. 

We will employ the Floquet formalism^'^ which has 
been extensively applied to study time dependent peri- 
odic evolutions in systems ranging from two level sys- 
tems, including simplified models of fiux qubits,— to 
more realistic molecular and nanoscaled systems^^. The 
Floquet method allows to transform the periodically time 
dependent Schrodinger equation into an equivalent in- 
finitely dimensional eigenvalue problem for a time inde- 
pendent Floquet matrix. In general several truncations 
schemes are employed in order to tackle the analytical 
solution and reduce the infinite Floquet matrix to an ef- 
fective finite dimensional matrix.— 

The paper is organized as follows. In Sec. |ll] we in- 
troduce the model Hamiltonian and equations for the 
Josephson flux qubit. In Sec. Hill we present numerical 
results for the amplitude spectroscopy for the Josephson 
flux qubit by direct numerical calculation and using the 
Floquet formulation for the time dependent Schrodinger 
equation in the case of an harmonic drive. In this section 
we compare our numerical results with recent experimen- 
tal realizations. In Sec. IIVI we propose a new amplitude 
spectroscopy method by changing the initial conditions. 
Numerical calculation using the Floquet formulation are 
presented. Finally, Sec. |V] contains a summary and a 
discussion of the most relevant points of our findings. 



II. MODEL FOR THE DEVICE FOR THE 
JOSEPHSON FLUX QUBIT 

The DJFQ consists on a superconducting ring with 
three Josephson junctions^, enclosing a magnetic fiux $ = 
/$o, with $0 = /i/2e, see Fig[TJ 




FIG. 1: Circuit for the DJFQ as described in the text. 
Josepshon junctions 1 and 2 have Josepshon energy Ej and 
capacitance C, and junction 3 has Josepshon energy and ca- 
pacitance Q times smaller. The arrows indicate the sign con- 
vention for defining the gauge invariant phase differences. The 
circuit encloses a magnetic flux <I> = /$o. 

The junctions have gauge invariant phase differences 
defined as (p\ , Lp2 and (^3 , respectively, with the sign con- 
vention corresponding to the directions indicated by the 
arrows in Fig[TJ Typically the circuit inductance can be 
neglected and the phase difference of the third junction 
is: (/93 = —(^1 + (^2 — 27r/. Therefore the system can be 
described with two dynamical variables: 1,(^2- The cir- 
cuits that are used for the DJFQ have two of the junctions 
with the same coupling energy, Ej,\ — Ej^i — Ej, and 
capacitance, Ci ~ C2 ~ C, while the third junction has 
smaller coupling £^,7.3 — aEj and capacitance C3 = aC, 
with 0.5 < a < 1. In terms of the two-dimensional coor- 
dinate If ~ {(fii, If 2), the hamiltonian of the DJFQ isi^ 

H = ^^lxa-'^^ + V{0), (1) 

where we have normalized % by the Josephson coupling 
energy Ej, and = %Ec/Ej, with Ec = e^/2C. The 
kinetic term of the hamiltonian corresponds to the elec- 
trostatic energy of the system, where the momentum op- 
erator is P = —ihV^p, with = (5^, 5^): and the 
"mass" tensor is given by the matrix m, 

_ / 1 + a -a \ 
I —a 1 -\- a J ' 

The potential term of the hamiltonian corresponds to 
the Josephson energy of the junctions, and is given by 

y(</?) = 2-t-Q;— cosiySi— cos(^2— acos(27r/+(y9i— 1^92) . (2) 
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Typical flux qubit experiments have values of a in the 
range 0.6 — 0.9 and 77 in the range 0.1 — 0.6i^"— 

In quantum computation implementations^ii^ii the 
Josephson flux qubit is operated at (static) magnetic 
fields near the half-flux quantum, / — 1/2 + 6f, with 
6f < 1. For a > 1/2, the potential of Eq.Q has the 
shape of a double- well with two minima (within the do- 
main — TT < ipi < TT, —77 < ip2 < tt). After a change 
of variables one can define the transverse phase ipt = 
{(fii +ip2)/2 and the longitudinal phase tpi = {ipi — <i52)/2, 
obtaining 

V{(fi) ~ 2 + a — 2 cos (fit cos (pi — a cos{2nf + 2(pi ) . (3) 

The two minima are along the longitudinal direction 
(fii, at {(pi,ipt) = (±(y5*,0) separated by a maximum at 
{ipi,ipt) = (0,0). Each minima corresponds to macro- 
scopic persistent currents of opposite sign. Experimental 
measurements are sensitive to the sign of the persistent 
current jiS. and therefore they detect the probability of be- 
ing on one side or the other of the double-well potential. 
For Sf > (Sf < 0) a ground state |— ) (|+)) with nega- 
tive (positive) persistent current is favored, with energy 
e± (X ±6f. At Sf = the two minima have the same 
energy, and the two lowest energy eigenstates (|5'o) and 
jvl/i)) are symmetric and antisymmetric superpositions of 
the two states (|— ) and |-l-))) corresponding to the macro- 
scopic persistent currents. To describe the dynamics of 
the device as a quantum bit, a two-level truncation of the 
Hilbert space is performed.^ In the subspace expanded by 
I^I^o) and l^*!), the hamiltonian of Eq. ([T]) is reduced to 

'Htls = ' (4) 

where Htls is written in the basis defined by |+) = 
(|*o) + |*i))/V2 and I-) = (|*o) - |*i))/V2. Here 
A = El — Eq is the two- level splitting at 5f = 0, 
and e = AnaEjSoiSf (for Sf < 1), with 6*01 = 
-(^-ol sin(2v3;)l*i) = -(+1 sm{2ipi)\+). (For typical val- 
ues of a and r/, one has 5*01 ~ 0.8 — 0.9). Most ex- 
periments control the system varying the magnetic field 
detuning Sf. The magnitude of the gap A depends ex- 
ponentially on a and rj. Recently it has been shown ex- 
perimentally that it is possible to manipulate the value 
of A by controlling a, replacing the third junction by an 
additional SQUID loop^^^s 

Landau-Zener-Stiickelberg (LZS) interferometry is 
performed by applying an harmonic field on top of the 
static field such that / f{t) with 

fit)^ f0 + fpSiTi{LUt) . (5) 

Hence Eq.(j3]) acquires an explicit dependence on time 
through e — > e{t) — eo+Asui{ujt) with eg — AnaEjSoiSf, 
(5/ = /q — 1/2 and A = AnaEjSoifp- The initial state 
corresponds to prepare the system in the ground state 
|0, /o) for the static field /q. 

For values of |(5/| ^ 1 and small driving amplitudes 
/p ^ 1, the DJFQ is adequately described as a two 



level system (TLS), whose time evolution under an har- 
monic drive does not have, in general, an exact solu- 
tion. Thus the dynamics is usually approximate d^^i^^ 
by free evolutions of the basis states mediated by non 
adiabatic Landau-Zener (LZ) transitions^'', with proba- 
bility Plz = exp (—277(5) with S = A'^/fpU. In the last 
case, explicit expressions for the occupation probability 
have been obtained in the fast (slow) driving regimCfiii^ 
(5< (>)1. 

Here we will consider the case of nearly fast driving, 
which corresponds to the series of experiments on ampli- 
tude spectroscopy performed in RefflSl For fo < 1/2, the 
system is started in the ground state |0, fa) « |-|-), which 
has a positive persistent current. In this type of experi- 
ments, one asks for the probability of switching to a state 
of negative persistent current: P|_)_^^|_^ = P-(t) during 
the time the harmonic pulse is applied. For the TLS in 
the fast driving regime, the occupation probability P_ is 
approximately given asJ^i^ 

P^^^it) = ^l|((l-cos(f}„0) , (6) 

r„ = A Jn{A/uj) , 

iln = \/ {nuj - eo)^ + r?i , 

being Jn{z) the order n Bessel function of the first kind. 
The resonance condition eg = is attained when the 
total phase accumulated over a single period of the driv- 
ing, Q — 27reo/w, satisfies Q = 27m for a given inte- 
ger riji^iii Under resonance, the occupation probability 
Pl^^it) 1/2(1 - cos(f7„t)) with = AJ„(A/a;). 
Notice that 0„ depends on the driving amplitude fp 
through A. 

Besides the time dependence, the average occupation 
probability is the key quantity for the spectroscopic anal- 
ysis performed in recent experimentsJ^ii^ In the case of 
a TLS in the fast driving regime, the average occupation 
probability obtained from Eq.® is a sum of Lorentzian- 
shape n~ photon resonances^ 



1 r2 

pTFD ^ ± iii /yN 

2^(nc.-6o)2+r2 • ^> 

Thus as eo (or Sf) is changed, different n-resonances are 
explored. In addition the Bessel function entering in r„ 
gives a quasiperiodic character to the patterns of res- 
onances as the amplitude fp is varied keeping the fre- 
quency uj flxed. 

The analysis of the positions of the resonances as a 
function of fp and Sf was the route followed in Refsimi 
in an effort to obtain the parameters characterizing the 
different avoided crossings of the flux qubit. 
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III. AMPLITUDE SPECTROSCOPY FOR 
COHERENT SYSTEMS 

A. Direct numerical calculation 

In this section we will focus on the study of the quan- 
tum dynamics of the D JFQ driven by the time dependent 
flux f{t) given in Eq.© for fo near 1/2 and varying the 
amplitude fp of the harmonic drive. 

In the absence of driving, i.e for fp ~ 0, the eigenvec- 
tors 5'„((^) and eigenenergies i5„ are obtained by solving, 



^n{0) = . (8) 



In Fig. ([2]) we plot the seven lower energy levels as a 
function of flux detuning /o, obtained by numerical di- 
agonalization of Eq. ^ using a discretization grid of 
Aip ~ 2tt/M and 27r-periodic boundary conditions on 
(f = {(pi,ip2). In this case we set rj — 0.25 and a = 0.8, 
close to the experimental values employed in flux qubits 
experiments.— The energy spectrum is rather sensi- 
tive to the values of rj and a, and in particular for 
the selected values, the energy landscape is quite in- 
volved, presenting many avoided crossings Aij in the 
range 0.45 < /o < 0.55. The slope of the energy lev- 
els dEn/dfo is proportional to the average current in the 
loop. Therefore an eigenstate with positive or negative 
slope, corresponds to a wave function mostly weighted 
in one side or the other of the double-well. A gap Aij 
opens at the avoided crossings of energy levels of oppo- 
site slope. We label the gaps Aij as the avoided crossing 
of the «-th level of positive slope with the j-th level of 
negative slope, see Fig.(|2l. ^This convention is different 
from the one used in Refll5l where a distinction among 
longitudinal and transverse modes is made in the labeling 
of the gaps.) 

It is evident from the energy level diagram of Figl2] 
that the description of the time evolution of the DJFQ 
in terms of a TLS is valid only for a very small range of 
amplitudes fp. 

In the presence of a finite driving amplitude fp, our 
first approach to the problem is to solve numerically the 
time dependent Schrodinger equation (we have normal- 
ized time by tj — h/Ej) 



' dt 



(9) 



We integrate numerically Eq. ([9]) with a second order 
split-operator algorithm, using a discretization grid of 
A(p = 27t/M and At — O.ltj. We use 27r-periodic bound- 
ary conditions on (p = [ipi, ip2). The system is started in 
the ground state |0, fo) for a given static field fo, obtained 
from the numerical solution of Eq. ([5]) . Experimentally, 
what is measured is the probability of being in a given 
state of positive, P+, or negative, P_, persistent current, 




FIG. 2: Lowest seven energy levels of the DJFQ as a function 
of flux fo for r] = 0.25 and a = 0.80. Arrows indicate the 
position of the avoided level crossings A^j measured from /o — 
0.497 (indicated by the vertical dashed line). Calculations 
were done using M = 256 (see the text for details). Energy 
is measured in units of Ej and flux in units of $o- 



which can be obtained as: 



P+it) = 1 - P.{t) = [ |^(^l,^2, 



t)\'^dipidip2 , 



(10) 
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FIG. 3: (Color online) Large amplitude spectroscopic dia- 
monds obtained for the Josephson flux qubit with rj — 0.25 
and a = 0.8. The intensity of P_ is plotted as a function of 
flux detuning 5f and rf amplitude fp. The system is driven 
at a frequency lu — 0.001 (in units of Ej/h). Numerical 
calculations were done by direct numerical integration using 
M = 128 and At — O.ltj. Data points correspond to a grid 
of A/o = 1 X 10^^ and Afp = 1 x 10"". 

where the integration is on one side of the double-well 
potential defined hy ipi > (i.e., — tt < t/Si < tt and 
'^2 ^ fi)- The quantity measured experimentally is 
the long time occupation probability, which is equiva- 
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lent to the time average probability in the stationary 
state. In Figl3] we plot the time averaged probability 
P_ = 1 — P+ as a function of the static flux /o and the 
amplitude fp of the harmonic excitation. The average 
is performed over several periods of the harmonic drive 
(typically ^ 20 — 100 periods, until convergence of the 
average). The plot is obtained by calculating points with 
a grid of A/o = 1 x lO'^ and Afp = 1 x lO"*. A pat- 
tern of "spectroscopic diamonds" is observed, similar to 
the one obtained in the experiments, which can be re- 
lated to the energy level spectrum of Figl5]as follows. At 
a fixed flux detuning Sf = fo — l/2, the first diamond, 
Dl, starts when the Aqo avoided crossing is reached, at 
fp = fis = /oo - /o = -Sf, with /oo = 1/2 the location 
of Aqo. The first diamond ends when the Aio crossing is 
reached at fp = fie = fo — fio, with /lo the location of 
AiQ. Then the second diamond, D2, starts when the Aqi 
avoided crossing is reached, at fp — f2s — foi — fo, with 
/oi the location of Aqi, etc. The spectroscopic diamonds 
of Fig O have much less contrast than in the experiments 
of Ref.lS, in which the contrast is due to fast intra-well 
relaxation. This induces population inversion, reducing 
the population in the sectors between the diamonds (i.e. 
for example between Dl and D2). In DJFQs with less 
relaxation effects, the picture should be closer to the one 
shown in FigI21 Within the first diamond a regular pat- 
tern of resonances can be qualitatively observed in FigO 
However to accurately describe all the resonances of Dl, 
as well as the complex structure of D2, a finer grid sam- 
pling A/o,A/p is needed. Furthermore, simulations at 
high fp, in the region of the second diamond and above 
(where the dynamics has more weight in higher energy 
levels), need a better discretization of the Schrodinger 
equation and averaging of the population for larger times. 
Therefore a finer description of the structure of diamonds 
needs large time consuming simulations of the full time 
dependent Schrodinger equation. Instead, in the follow- 
ing we will employ an approach based on the Floquet 
formalism, more adequate for time periodic hamiltoni- 
ans. 



B. The Floquet Formulation 

For a finite driving f{t), we write V. = Ho + SV{0, t) 
with Ho corresponding to Eq. ([T]) with f = fo (i.e. the 
time independent part of the hamiltonian) and 

6V{(p, t) = aEj sm[27r facit)] sin(27r/o + ifi - 1P2) + 

2aEj sm^[nfac{t)] cos(27r/o + ipi - (^2) (H) 

where fac{t) — /psin(aji). The potential defined in 
Eg. (1111) is periodic in time with period T = 27r/w. Thus, 
according to the Floquet theorem,-- the time dependent 
Schrodinger equation ^ has a solution that can be writ- 
ten as 



where ^a{t) = ^a{t + T) and Sa is known as the 
quasienergy or Floquet eigenvalue. 

Substituting expression into Eq.® we obtain an 
eigenvalue equation for the quasienergies: 



HF{t)<i>o.{t) = eMt) , (13) 
where the Floquet Hamiltonian is defined as 



Hp(t)=H{t)-i 



d_ 

di 



(14) 



Since the function $a(0 i*^ periodic in time it can be 
expanded in a Fourier series. We introduce the standard 
Floquet nomenclature's, and define \n,k) — \n) (g) \k), 
where n is an index that labels the eigenstates of Ho and 
fc is a Fourier index. Then 



iktot 



(15) 



k — — OQ 



where {n,k\4)a) is a Fourier amplitude. From Eq. ()13lll5|) 
it is straightforward to write 



ea{n,q\4'a) = '^'^{n,q\HF\'m,k){m,k\(j)a) (16) 

m k 

where Hp is the Floquet Hamiltonian previously defined 
whose matrix elements are given by 



(n, q\HF\m, k) = (En + quj)5m,n5k,q + 



2ti 



{n\W{ip,t)\m)e-'^'i-^'^'^'dt 



(17) 



Then the time dependent problem is reduced to solve the 
eigenvalue equation Eq. (fT6|) . 

We need to calculate the matrix elements of the Flo- 
quet Hamiltonian defined in Ea. (|T7l) . As it was men- 
tioned before, the first term in this equation is obtained 
by numerical diagonalization of Eq. (jS]), using a dis- 
cretization grid of A(^ — 2t:/M and 27r-periodic bound- 
ary conditions on ~ ((^1, Lp2). The second term in Eq. 
(fTT]) is written as 



Jo 



2tt/u 



(18) 



where SV{t) was defined in Eq. pTj) and I — q — k is a.n 
integer. The integration is straightforward and we obtain 



*a(i) 



(12) 
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where 5„„ = (n| sin (27r/o + v^i - </?2)|?ti), Cnm = 
{n\ cos (27r/o + — (p2)\iTi) and Ji{x) is the Bessel func- 
tion of first kind of order /. 

Then we have all the ingredients to construct the Flo- 
quet matrix 

{n, q\HF\m, k) = {En + qLo)5,n,n5k,q + KL'' , (20) 

where q and k range over all integers form — cx) to oo. In 
order to solve the problem numerically we must truncate 
the Floquet matrix, Eg. (1^01) . The truncated matrix is 
of dimension Nd = {2K + l)Ni where K is defined by 
the maximum value of the Fourier index and Ni by the 
number of levels considered in the diagonalization of Eq. 

Floquet eigenstates and quasienergies contain all the 
information to construct the large amplitude spectro- 
scopic diamonds. Following the experiments, we take as 
initial state the ground state of 'Ho for a given value of 
flux detuning /q that here for simplicity, we denote |0). 
The initial state is prepared at a time io? and then at a 
time t the evolved solution |5'(t, io)) can be expanded in 
the basis of eigenstates of Ho as 

■^{ipi,Lp2,t,t(i) ^ ^C„(<,to)Xn(<y3l,'P2) (21) 
n 

where Xn{'~Pi,'^2) = {'^iT'f2\'n) is the wave function rep- 
resentation of eigenket |n) in terms of the variables 

(<Pl,'<f2). 

Using the Fourier amplitudes of the Floquet eigen- 
states {n,k\4>a), and their corresponding eigenenergies 
£a, one obtains the coefficients^^ 

Cn(i,to) = EE("'^l'^^)<<^^|0'0)e""'^*"'''^e^"'* (22) 

k /3 

which are the probablity amplitudes that the system ini- 
tialy in the ground state at time io evolves to a state \n) 
by time t according to the time-periodic Hamiltonian. 
This equation can be interpreted as the amplitude prob- 
ability that the system initially in the Floquet state |0, 0) 
at time evolve to the Floquet state |n, k) by time t ac- 
cording to the time independent Floquet Hamiltonian, 
sumed over k with weighting factors exp (icukt) . 

We can now calculate the time dependence of the prob- 
ability P+ (or P_ = 1 — P_(.) of a measurement of a 
positive (negative) state of persistent current, replacing 
Eqs.(I2Il) and (HH) into Eq.dfOl), and obtaining, 

p+{t,to) = y^y^ A„„t(t,^o)p„„ , (23) 

n rn 

where the coefficients 



Pnm = / Xni'Pl,'f2)Xm{'^l,f2)d'fldip2, (24) 

Jw 



are evaluated by numerical integration using the eigen- 
states of Eq.® and W is the triangular sector of the 
two dimensional space defined hy ipi > (i.e., — tt < 
(ySi < TT and tp2 < fi)- The time dependent coefficients 
Kim(t,to) = c„(t, to)c*„(i, io) are calculated as 

A„m(t,io) = Efc,(E/3,'y('^>^l<^/3)(</'^i|0>0) 

(O,O|0^)(0^|TO,Oe-*'^^-^'^)(*-*o)e*'^(''-''* . (25) 

We can calculate now the time average of P+ (t) . To this 
end, we average Xnm(t,tQ) over many periods consider- 
ing that for long times e-tl'-'ik-i)+ey-e^]t _ S^i^k_i)^i;^_^^. 
Using the periodic properties of the quasienergies and 
eigenvectors of the Floquet Hamiltonian^ we get 

Km{to) = EfcjE/3("'fc|'/'/3)(</'/3|0:0) 

(0,fc- ?|(?!)0)(</./3|m,fc)e*"('=-')*« (26) 

In experiments the initial time, or equivalently the initial 
phase of the field seen by the system in repeated realiza- 
tions of the measurement, is not well defined. Then, the 
quantity of interest is the transition probability averaged 
over initial times,— in which case e*"('^~')*o = Sk,i- Fi- 
nally, we obtain 

P+ = y^PwrnAnm , (27) 
n,m 

where 

Anm = I] ('/'/3 |0, 0) (0, 0\^p) {cjjp |m, k) . 

fi k 

Thus once equiped with the Floquet quasienergies and 
eigenstates one can either compute the time dependent 
occupation probability Eq. (j23p or the time average prob- 
ability Eq. (f27| ■ which is indeed an exact average. 

To summarize, we follow this procedure: (i) We take 
the parameters 77 and a from the experimental device. 
(Here we use 77 = 0.25 and a = 0.8.) (ii) For each value 
of the magnetic flux /o wc solve numerically the eigen- 
value equation ([S]) obtaining the eigenstates Xn{vi,f2) 
and eigenvalues i?„, with a discretization A<y9 = 2tt/M 
(M = 256 - 1024). (in) We evaluate the Floquet matrix 
elements of Eq. ([TO)) and the coefficients Pnm of Eq. . 
(iv) We solve numerically the Floquet eigenvalue equa- 
tion, Eqs. ((T6|) and (|20l) obtaining the Floquet eigenvalues 
and the components (n, gl^^) of the corresponding 
eigenvectors. A truncation of the matrix is performed: 
we consider Ni energy levels (results for different Ni will 
be shown) , and we consider 2K + 1 Fourier components 
m —K <k < K. Large K is chosen until convergence of 
the quantities of interests. (X ~ 150 — 400) (v) The time 
dependent occupation probability Eq. (P5)) for <o = 0, or 
the time average probability Eg. ([27)) are then evaluated, 
with the sums over energy levels between and Ni — 1, 
the sums over k components between —K and K , and the 
sums over Floquet states between 1 and Nd = Ni{2K+\). 
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1. Time Dependent Occupation Probabilities 
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FIG. 4: (Color online) P_ as a function of time r — ijjt/2n for 
77 = 0.25 and a = 0.8. (a) eo = t^, for a; = 0.001, 5/ = -0.11 x 
10"^, and fp = 0.24 x 10"^ (black solid line). The fast driving 
TLS Eq.© for n = 1 with f^i = Aoo Ji(47rS'oi/p/w) = 1.87 x 
10~* is plotted for comparison (red dashed line). (b)eo ~ 3a; 
for Lj = 0.001, Sf = -0.33 x 10"^ and fp = 0.35 x 10"^ (black 
solid line). The fast driving TLS Eq.® for n = 3 with fis = 
Ago J^3(47rS'oi/p/a;) = 1.14 x 10~* is plotted for comparison 
(red dashed line), (c) eo = 3cj for a higher frequency lj — 
0.003, 5f = -0.99 X 10"^ and fp = 1.05 x 10"=^ (black solid 
line). Notice that the oscillations on time scales --^ tt/cj are 
smaller for u = 0.003, as described in the text. The fast 
driving TLS Eq.(|6]) for n = 3 is plotted for comparison (red 
dashed line), (d) Out of (n = 3) resonance, for cj — 0.001, 
5f = -0.32 X 10"^. fp = 0.35 x 10"^ (black solid line). Notice 
that the oscillations on time scales ~ tt/oj persist. Numerical 
calculations were done with A'^; = 6 levels, K — 180 and 
M = 1024. Calulations with Ni — 2 levels overlap almost 
exactly with the Ni = 6 calculations in this case. 

We start analyzing the explicit time dependence of the 
probability P-{t) after a driving of duration t is applied. 
To this end, we evaluate numerically Eg. ([25)) for to = 0. 
Different initial states shall correspond to prepare the 
system in the ground state \n — 0; /o) for different /o < 
0.5. 

For values of (5/ = /o — 1/2 and (rather small) driv- 
ing amplitudes fp, such that Ago is the only relevant 
avoided crossing (see FigH]), the DJFQ can be described 
as a TLS. In Figglwe plot P_(t) 1 - P+(r) as a func- 
tion of time T = u)t I 2t: in units of the pulse period, for 
small fp. The numerical calculations were performed 
with the Floquet formalism, employing in Eq. (P5| the 
lowest six energy levels [Ni — 6) and K — 150 — 250. 
For this case, we find that calculations with Ni — 2 lev- 
els overlap almost exactly with the Ni = 6 calculations, 
and can not be distinguished in the plot, since the TLS 
approximation is correct for small fp, as expeceted. In 
panel (a) we consider the case with frequency uj = 0.001 



and with Sf = —0.11 x 10~^ such that it corresponds 
to the fast driving TLS n — 1 resonance, eg — ^- For 
the selected amplitude, fp — 0.25 x 10^'^, the behav- 
ior of P- (r) is on the global scale rather well described 
by the fast driving approximation for a TLS given by 
Eq. ([S]) . The numerically obtained frequency is very close 
to ill ~ Ti = 1.87 X 10~*, in agreement with the on - 
resonance relation written in Eq.(|n]). However, P-ij) ex- 
hibits sudden jumps mediated by additional oscillations 
with n local maxima on time scales ^ tt /lo. These oscil- 
lations reflect the quantum mechanical interference be- 
tween consecutive passages through the avoided crossing 
AoQ. This behavior is not captured by the fast driving 
TLS expression Eq.([6|) plotted for comparison by the red 
dashed line. In the present case 5 — Ago/w fp ~ 0.5 and 
then, the interference effects are important. In panel (b) 
we consider for the same frequency as in (a), the case 
of 6f = -0.33 X 10"^ to select the fast driving TLS 
n = 3 resonance (eg = 3a;). For a small amplitude 
fp — 0.35 X 10^'^, the qualitative description given for 
panel (a) holds. For the driving parameters used in panel 
(b) we have an adiabaticity parameter i5 ~ 0.25. In this 
case the oscillations on time scales ~ tt/w, exhibit n = 3 
local maxima. In panel (c) we also consider the case of 
the n — i resonance but at a higher frequency, uj — 0.003, 
for which 5 ~ 0.03. In this case the fast driving approx- 
imation is more adequate, the local oscillations on time 
scales ^ tt/w are washed out, improving the agreement 
with the fast driving TLS expression Eq.(l6l). This can 
be understood taking into account that the (adiabatic) 
Landau- Zener transition probability^! at a single avoided 
crossing diminishes as the frequency oj increases. In Fig|4] 
(d) P_ (r) is plotted for 5f = -0.32 x lO"^ and w = 0.001, 
a value which is out but close to the n — i resonance. 
We employ the same amplitude fp = 0.35 x 10^"^ as in 
FigHl (b). Notice that in this case is max[P_(r)] < 1, 
as expected in the off-resonance situation, but the short 
time scale oscillations on time scales ~ tt/w persist. 

In Figl5]we show that for a larger amplitude /p = 16 x 
10~^ the TLS approach breaks down. For this amplitude, 
the system is driven close to the avoided crossing Agi 
(see Figl2) . In panel (a) we show the case where eg = 3a; 
is satisfied while in panel (b) we show the case out of 
the eg = 3a; condition. As expected, the magnitude of 
max[P_ (r)] in this case is completely unrelated with the 
resonance conditions observed at smaller fp. We also 
compare the results obtained considering up to Ni — Q 
levels when evaluating Eq. (P^ with the case with only 
Ni = 2 levels. As it is evident in the plots, for this large 
amplitude more than two levels are needed to describe 
the behavior of P_ (r) , since most of the population is at 
the higher energy eigenstates. 



2. Average Occupation Probabilities 

Besides the time dependence, the average occupation 
probability Eq. p7)l is the key quantity for the spectro- 
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FIG. 5: (color online) P_ as a function of time r = ujt/2-n for 
rj — 0.25 and a — 0.8. The flux qubit is driven at a strong 
amplitude with /p = 16 x 10~^. (a) resonance condition for 
the TLS in the fast driving regime: eo = 3a; for uj = 0.001, 
Sf = -0.33 X 10"^ (b) eo / mo, for to = 0.001, Sf = 
—0.32 X 10"''. Numerical calculations were done with TV; = 6 
levels, K = 300 and M = 1024. Calculations with Ni ^ 2 
levels are plotted for comparison (red line). 



scopic analysis performed in recent experimentsJ^ii^ We 
analyze the patterns of P_ in the {fp,6f) space. As we 
shall show below, as the amplitude fp is increased, P_ 
will exhibit a reacher and more involved structure. 

Regarding the parameters employed in the numerical 
calculations, for small driving amplitude fp and flux de- 
tuning /o ~ 0.5, the discretization grid needed to com- 
pute the eigenstates |n; /o) can be constructed with quite 
small values of M ~ 1281^^ However, in the strong driv- 
ing regime, we need to use M = 256 — 1024. In addition, 
we employ K ^ 150 — 400 to attain convergence in the 
values of P_ . As we already mentioned TV; is mainly de- 
termined by fp. All the calculations, otherwise specified, 
have been performed with the lowest six levels. 

In Fig. mja) we show the contour plot of P_ as a 
function of {fp,Sf) for J/ < and a; 0.001. The 
range of values of fp and Sf has been selected to explore 
the region of the energy spectrum of Figj^] containing 
the avoided crossings Aqq, Aiq, Ah and A20. The plot 
is obtained by calculating points with a grid of A/o — 
2 X 10"^ and Afp = 3 x 10~^. A clear pattern of maxima 
and minima forming a half diamond-like structure can 
be observed in Fig. IUJa). For the sake of clarity we 
have drawn lines indicating the boundaries of the first 
(half) diamond Dl and the beginning of the second one 
D2. For Sf > the diamond pattern is completed with 
1 — P_. The qualitative agreement with the diamond 
structure observed in Ref. |l5|is evident. However, in this 
coherent regime the spectroscopic diamonds have much 
less contrast than in the experiments of Reffisl 

In Fig. ini[b) we show a similar contour plot of P_ 
for a higher frequency uj = 0.002. A very similar pat- 
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FIG. 6: (color online) Large amplitude spectroscopic (half) 
diamonds obtained for the Josephson flux qubit with rj — 0.25 
and a = 0.8. The intensity of P- is plotted as a function 
of flux detuning 5f and rf amplitude fp. (a) The system is 
driven at a frequency lj — 0.001. Calculations were done using 
TV; = 6 (six levels), 150 < K < 300 and M = 1024. The black 
solid lines indicate the edges of the flrst diamond Dl and a 
the beginning of the second one D2. The black dot indicates a 
particular value of flux detuning 5f ~ —0.003 for fp = 0. See 
text for a detailed analysis, (b) The flux qubit is driven at a 
frequency uj = 0.002. Calculations were done using TV; — 6, 
150 < K < 250 and M = 1024. Data points correspond to a 
fine grid of A/q = 2 x 10"^ and Afp = 3 x 10~^ 



tern of maxima and minima forming a half diamond-like 
structure can be observed in Fig. (HKb). For this larger 
driving frequency the distance between resonances is in- 
creased, as it is expected from Eq.([6|). It is important 
to remark that when we increase the driving frequency 
and thus the sweep rate, we lose resolution in the ob- 
tained spectroscopic diamonds as it can be checked by 
inspection of FiglHlJa) and FigEl^b). 

Emulating the experimental protocol, in the following 
we analyze the structure of the diamonds in order to ex- 
tract spectroscopic information, focusing on a fixed fre- 
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quency, ui = 0.001. We start by considering a particular 
static flux detuning /o ~ 0.49691 i.e. df ~ -0.00309 
(black dot in FiglBJa) and vertical dashed-dotted line in 
FigHI) that satifies for uj = 0.001 the fast driving TLS 
n — 28 resonance condition, eo = 28a;. The initial ground 
state is |0; 0.49691) and thus for /p = is P_ = 0. 
As the driving amplitude is increased, the net transfer 
of population over many driving periods will translate 
in a finite value of P_. Along the horizontal line de- 
fined at 6f ~ -0.00309 in Fig. [DJa), the first diamond 
Dl starts at ~ 0.003 {D labels the parameters ex- 
tracted from the diamonds). From Fig 15] one can check 
that this value is roughly the threshold amplitude needed 
to reach the first avoided crossing Aqo for the considered 
value of /o = 0.49691. For fp > 0.003 the muhiple pas- 
sages through the avoided crossing Aqo are reflected in 
the observed interference pattern that, up to fp ^ 0.01, 
it is rather well described by the quasiperiodic behavior 
of the Bessel function Jn=2a{4,TTaEjSQifp/u!) entering in 
the definition of r„=28 in Eq.Q. If the driving amplitude 
fp is further increased the interference patterns persist, 
but the positions of the maxima and minima do not fol- 
low the Bessel function dependence, as the description of 
the resonances in terms of the fast driving TLS formula is 
not accurate for these amplitudes. A detailed evidence of 
this behavior is shown in FiglTl^a) where we plot a cut of 
P_ along the horizontal line Sf = —0.00309 depicted 
in Figl6ja), corresponding to the condition £o — 28aj 
(n = 28 resonance condition for the TLS in the fast driv- 
ing regime), together with the same quantity computed 

(2) 

keeping only the lowest 2 levels, P_ , i.e. with Ni = 2. 
In FiglTl^b) we plot a case slightly different (slightly-off 
n = 28 resonance condition for the TLS in the fast driv- 
ing regime), for Sf = —0.00308, where P_ is clearly 
smaller. For small amplitudes, P_ is nicely followed by 

(2) 

P_ , that indeed reproduces quite accurately the inter- 
ference patterns both on-resonance and off-resonance, in 
FigEfa) and (b), respectively. However, for fp > 0.01 the 

(2) 

departure of P_ from the actual behavior is notorious, 
even before the end of the first diamond Dl. In addition, 

(2) 

P_ — ^ refiecting the fact that higher levels besides the 
lowest two are populated as fp increases. 

Unlike its beginning, the end of Dl is not so sharply 
defined, showing a rather poor change in contrast. In the 
experiments of Ref .T^ the data show a larger reduction in 
contrast, due to fast relaxation through intra- well transi- 
tions. Furthermore, in the experiment is Aio S> Ago and 
the population transfer is dominated by the transition 
at the avoided crossing Aio, with no explicit signatures 
of additional multiple passages through the extra avoided 
crossing An (see Fig. 1(c) in Reffish. In our case we have 
verified that, affhough Aio = 2 x 10"^ > Am = 3 x 10"" 
is roughly the same relation than in the RefllSl additional 
transitions at An = 1 x 10~^ contribute to sustain the 
values of P_, resulting in an effective reduction of the 
contrast at the end of Dl. It is plausible that in the ex- 
perimental qubit^^, a high value of An gives, unlike our 




FIG. 7: (color online) P_ (black solid line) calculated employ- 

ing the lowest 6 levels and P_ (red dashed line) calculated 
employing the lowest 2 levels as a function of the driving am- 
plitude. Calculations were done using K — 150 — 300 and 
M = 1024. (a) 5f = -0.00309 and (b) 5f = -0.00308. The 
circle, square and triangle denote the beginning of the first 
diamond Dl, the end of Dl, and the beginning of the second 
diamond D2 respectively, for the present value of 5f. See text 
for details. 



case, a negligible transition probability at this avoided 
crossing. 

Here, we find that despite the poor contrast in a highly 
coherent DJFQ, the end of the first diamond is quite iden- 
tifiable at /j^ - 0.013, giving a value of //^ = 0.484 in 
good agreement with the position of the second avoided 
crossing Aio obtained from the analysis of the spectrum 
depicted in FiglJJ Thus, the boundaries of the first di- 
amond give a rather satisfactory determination of the 
position of Aoo and Aio, respectively. 

For values of fp > 0.013 the competition between the 
different transitions at the avoided crossings Aoo, Aio 
and All, turns the interpretation of the pattern followed 
by P_ rather complicated. However the beginning of the 
second diamond D2 at f£. — 0.019, gives the position of 
the avoided crossing Aoi, at /J^ = 0.516 very close to the 
exact value (see Figl2|). 

In analogy with the experimental analysisiS,, the dia- 
monds profiles can be studied in more detail for a given 
amplitude fp and sweeping the flux detuning /o. To this 
end, we select two vertical lines in Fig|6]that correspond 
to fp = 0.001 and fp = 0.018, respectively. 

In the upper panel of FiglS] we show P_ along fp — 



0.001. A dashed line indicates the value /( 



Dl 



0.4988 



at which the vertical line defined by fp = 0.001 inter- 
sects the lower edge of the Dl (see FiglH]). As expected, 
the symmetry P_ — >• 1 — P_ around fo = 0.5 holds. 
The proflle of P_ , sweeping the different equally spaced 
n-resonances as fo changes, seems to be in very good 
agreement with the predicted TLS resonance pattern. 



10 




FIG. 8: (color online) P_ (black solid line) as a function 
of the static flux fo for two different values of the driving 
amplitude fp — 0.001 (upper panel) and fp = 0.018 (lower 
lanel) calculated using Ni = 6. In both panels, a dashed line 
denotes the value of fo that for fp = 0.001 {fp — 0.018) gives 
the beginning of the first(second) diamond Dl (D2). As a 

(2) 

comparison we plotted in red dashed line the P_ computed 
employing the lowest 2 levels. Numerical calculations were 
done with AT = 150 - 300 and M = 1024. See Fig. [6] and text 
for further details. 



— (2) 

Indeed we have included the results for P_ computed 
employing the lowest 2 levels, which are essentially su- 
perimpossed to the actual P_. An estimate of the ob- 
served number of resonances obtained employing Eq.([7]) 
is n = [eo/w] — [4:TTaEjSoi\Sf\/uj], being [...] the integer 
part. In this case is \Sf\ = |0.5 - f^^\ = 0.00122 and 
AwaEjSoi ~ 9. Thus for oj = 0.001 we obtain n = 11, 
which is exactly the number of maxima displayed in the 
upper panel of Fig l8] for the selected range of flux detun- 
ings. Thus this analysis, complemented with the previous 
one performed in Fig (71 confirms that close to the begin- 
ning of the first diamond Dl, the TLS description is quite 
accurate. 

The lower panel of FigE displays P_ for fp = 0.018 
and the beginning of the second diamond D2 is indicated 

/[P^ = 0.4974 by the vertical dashed line. The erratic 
pattern of resonances in P_ is in correspondence with 
the results presented in Fig 17] for amplitudes inside the 
second diamond D2. As it occurred in that case, the 
competition between the transitions at different avoided 
crossings gives a profile of the occupation probability that 
strongly departs from a simple interference pattern as 

given by Eq.([7l) and/or for the pattern displayed by P_ . 

Coming back to complete the spectroscopic analysis of 
FiglH the end of the second diamond should be expected 
at fp ~ 0.018 corresponding to the position of the A20 
avoided crossing. However, as it is easily checked from 
Figini we do not obtain the end of the second diamond for 
this value of fp. Indeed D2 starts for a larger value of the 
amplitude. In our case A20 ~ 1 x 10~^ ^ Aqq and there- 



fore the transition probability cx A2o^/w/p4 —^0. As a 
consequence, the spectroscopic diamond contains no vis- 
ible information on the avoided crossing A2o- This small 
gap should correspond to a crossing of transverse modes, 
which are not easily probed by the driving f(t) which 
acts mainly along the longitudinal ipi direction. Similar 
drawback for detecting transverse modes has been re- 
ported in the experiment A possible way to increase 
the resolution of this gap is to increase the transition 
probability by reducing the driving frequency ui. How- 
ever, larger driving periods could be concomitant with 
the loss in resolution of the indivual n-resonancesJ^. 

Obtaining numerically the diamond pattern besides 
the beginning of the second one D2 is a formidable task, 
essentially due to the extremely large CPU needed. Fig. 
ini shows a cut of P_ for Sf = —0.0020 up to ampli- 
tudes fp ^ 0.04. For the larger amplitudes {fp > 0.02), 
we needed to perform the calculations employing 8 levels 
{Ni = 8). For the small amplitudes the edges of the first 
diamond Dl, and the beginning of the second diamond 
D2 are clearly visible in the abrupt changes exhibited 
by P_. On the other hand, as we have already men- 
tioned, no evidence of population transfer is obtained for 
the transverse avoided crossings A20 and Ao2. For larger 
amplitudes, one can distinguish a region where P_ is 
small (P_ < 0.1) as separating the end of the second 
diamond D2 and the beginning of the third diamond D3. 
In Fig. ini we show the points in fp where the avoided 
crossings at A30 (triangle down) and A03 (open circle) 
are reached. 



P 0.3 




FIG. 9: (color onhne) P- for Sf — —0.002 calculated employ- 
ing the lowest 8 levels as a function of the driving amplitude. 
Calculations were done using K = 150 — 400 and M = 1024. 
The red symbols in the axis show the values of fp for which 
the different avoided crossings are reached (see Fig. 2). The 
vertical dotted lines indicate the borders of the different spec- 
troscopic diamonds Dl, D2, D3 
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IV. EXCITED STATE AMPLITUDE 
SPECTROSCOPY 

In this section we explore an alternative amplitude 
spectroscopic method to study the quantum dynamics 
of the DJFQ starting from a different initial condition. 
The gedanken experiment consist on preparing the sys- 
tem in the first excited state \n = 1, /o) for /o < 0.5. In 
this way, and depending on the value of /o chosen, the 
Aoo or Aio avoided crossing could be reached first as the 
amplitude is increased. With this initial state, the time 
averaged occupation probability is: 



n,m— 



(28) 



where 



K^^Y. E (",fc|'/'/3>('/'/3|l,O)(l,O|0^)(0^|m,fc) . 

p=l k=-K 

If we compare with Ea. (l?7|) we see that now de- 
pends on the amplitude (0^|1,O) instead of (0^|O, 0). In 
general, in highly coherent devices, one can define an 

average occupation probability P_|_ depending on the 
initial state |s), with the amplitude (0^|s,O) instead of 
((/)^|0,0) inEq.dlZl). _ 

For values of /o < 0.5 and /p we have P+ = 1 
(P_ = 0) and P*; = (pI = 1). In Fig. [TU] we 
plot P_ and 1 — P _ as a function of the driving ampli- 
tude for a fixed value for the detuning {5f = —0.00066), 
such that the avoided crossing Aqo is reached for smaller 
amplitudes than the Aio- For small values of fp the 
system behaves as a TLS and both probabilities give 
the same information. As the driving amplitude is in- 
creased {fp ~ 0.015) noticeable differences between P_ 
and 1 — P_ emerge. Indeed, in the highly coherent case 
the initial condition plays an important role in the quan- 
tum dynamics of the system already when approaching 
the second avoided crossing Aio ■ 

In Fig. [TT] (a) we show the contour plot of 1 — P_ as a 
function of (/p, 5f) for 5/ < and w = 0.001. The range 
of values of fp and 5f and the grid are the same used to 
obtain the results showed in Fig. [51 By inspection of Fig. 
ini (a) and [11] (a) we conclude that both amplitude spec- 
trospy methods give the same information before the Aio 
avoided crossing is reached. While the end of the first di- 
amond Dl is determined with a good contrast with the 
excited state amplitude spectroscopy method, the begin- 
ning of the second diamond D2 is still very difficult to 
determine. 

In order to analyze the information hidden in the dif- 
ferent espectroscopic diamonds we plot the difference be- 
tween the probabilities used to construct the diamonds 
in Fig. |n](a) and Fig. [n](a). In Fig. [II](b) we show the 
contour plot of P_ + P _ — 1 as a function of (/p, 5f) for 




FIG. 10: (color online) P_ (red line) and 1 — P1 (black line) 
as a function of the driving amplitude for 5f = —0.00066. 
Calculations were done employing the lowest 6 levels using 
/S' = 150 - 300 and M = 1024. 



6f < and lu — 0.001. As we mentioned, the difference 
is zero for values of fp such that the Aio is not reached. 
It is interesting to mention that now the end of the first 
diamond Dl and the beginning of D2 can be determined 
with a rather good contrast, due the cancellation of some 
intrincated interference patterns present in both P_ and 
1-Pl. 

Finally we evaluate the probability P_ as a function 
of the driving amplitude for a fixed value of the flux de- 
tuning close to the Aio avoided crossing. In this case P_ 
follows the TLS behavior up to the Aio avoided cross- 
ing, but once the driving amplitude reaches the Aoo ad- 
ditional levels should be included in order to properly 
described the quantum dynamics of the DJFQ. In Fig. 
[T^ we plot P_ for /o ~ /lo = 0.484, as a function of 
the driving amplitude up to values that drive the system 
close to the A30, for which we have to employ eight levels 
in the numerical calculations. In the figure the changes in 
P_ reveal the position of the different avoided crossings 
allowing a clear detection of Aio, Aoo and A30. 



V. SUMMARY AND CONCLUSIONS 

We have numerically solved the quantum dynamics of 
the device for the JFQ under strong harmonic driving 
in the fully coherent regime. Starting from the ground 
state we have studied the temporal evolution of the oc- 
cupation probability and analyzed the spectroscopic di- 
amonds obtained for the time average occupation as a 
function of flux detuning and driving amplitudes. We 
have shown that for small amplitudes the description in 
terms of a TLS reproduces very well the observed pattern 
of Landau-Zener-Stiickelberg interferences, as expected. 
On the other hand, the TLS description breaks down for 
driving amplitudes such that the avoided crossing Aio 
is reached. The spectroscopic diamonds exhibit in this 
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FIG. 11: (color online) Large amplitude excited state spectro- 
scopic (half) diamonds obtained for the DJFQ with r) = 0.25 
and a = 0.8 at a driven frequency uo = 0.001. Calculations 
were done employing TV; = 6 (six levels), 150 < < 300 and 
M = 1024. (a) Intensity plot of 1 - p1. (b) Intensity plot of 

p_-(i-r_) 
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FIG. 12: (color online) P1 calculated employing the lowest 8 
levels as a function of the driving amplitude for 5f = —0.014. 
The system is driven at a frequency uj = 0.001, K = 150 — 300 
and M = 1024. The red symbols show the position of the 
different avoided crossings. 

excited state amplitude spectroscopy, as performed in 
FiglTlTb). can now bring good contrast for the resolution 
of the second diamond. In a perfectly coherent closed sys- 
tem one could continue even further, performing another 
amplitude sweep starting from the second excited state, 
compare it with the results obtained starting from the 
first excited state, and so on. Of course, in a real system 
the possibility of these "excited state amplitude spec- 
troscopies" will be strongly limited by decoherence and 
relaxation processes. In current highly coherent DJFQ 
(with dcphasing times of the order of 1/xs) the first ex- 
cited state amplitude spectroscopy seems to be feasible. 
In this case, this could give an important indication of the 
coherence of the device as wells as additional and com- 
plementary information of the multilevel structure of the 
energy spectrum of DJFQ. 



case interference patterns with a rather complex struc- 
ture, due to the coherent evolution among all coupled 
energy levels. This situation is different from the exper- 
iment of Reffisl where there is a higher contrast in the 
diamond patterns due to the intra-well relaxation and 
short coherence times. In spite of this, in the fully co- 
herent regime explored in this work, we find that the 
edges of the diamonds clearly define the position of the 
different avoided crossings, as can be observed in FiglHl 
for example, even when the contrast is rather poor. In 
Section IIVI we have proposed a way to obtain further in- 
formation in this case. In a highly coherent DJFQ it is 
possible to prepare the system in the first excited state, 
for example with a tt Rabi pulse. From there, the ex- 
cited state amplitude spectroscopy could be performed. 
A comparison of the occupation probabilities obtained 
from the ground state amplitude spectroscopy and the 
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